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ASCA and Radio/RXTE Observations of GX 339-4 


M.A. Nowak, J. Wilms 1 , J.B. Dove 2 , R.P. Fender 3 * 
JILA , Campus Box 440, Boulder , CO 80309-0440, USA 

Abstract. We have analyzed three separate archival Advanced Satellite 
for Cosmology and Astrophysics (ASCA) observations and eight separate 
Ross : X-ray Timing Explorer (RXTE) observations of the black hole can- 
didate GX 339—4 in its low luminosity, spectrally hard state. Three of 
the RXTE observations were strictly simultaneous with 843 MHz and 8.3- 
9.1 GHz radio observations. All data sets show evidence for an « 6.4 keV 
Fe line with equivalent widths « 20-100 eV. ‘Reflection models’ show 
a hardening of the RXTE spectra with decreasing X-ray flux; however, 
these models do not exhibit evidence of a correlation between the photon 
index of the incident power law flux and the solid angle subtended by the 
reflector. None of the models fit to the X-ray data, however, simultane- 
ously explain the observed radio properties. We argue that the spatial 
extent of the observed radio emission is at least O(10 7 GM/c 2 ). Timing 
analysis reveals that all observations save one show evidence of a persis- 
tent /qpo « 0.3 Hz QPO. The broad band (10 _3 -10 2 Hz) power appears 
to be dominated by two independent processes that can be modeled as 
very broad Lorentzians with Q £ 1. Similar to Cyg X-l, the hard pho- 
ton variability is seen to lag the soft photon variability with the lag time 
increasing with decreasing Fourier frequency. The magnitude of this time 
lag is seen to be positively correlated with the flux of GX 339-4. 


1. Introduction 

GX 339-4 is a persistent, but highly variable, galactic black hole candidate that 
shows a wide variety of spectral states. It has an undetermined mass, and lies at a 
probable distance of « 4kpc (Zdziarski et al. 1998, and references therein). Here 
we discuss spectrally hard, low luminosity (3-9 keV flux £ 10” 9 ergs s” 1 cm” 2 ) 
states. We have analyzed three archival ASCA observations and a series of eight 
RXTE observations, three of which were simultaneous with radio observations 
(also discussed in Corbel et al. 1998 and Hannikainenen et al. 1998). Detailed 
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Figure 1. Reflection fraction vs. photon index, T, for models fit to 
PC A data only (squares), and models fit to PC A and HEXTE data 
simultaneously, but allowing the photon indeces and normalizations to 
be different (diamonds; HEXTE photon index shown). Middle: Photon 
index, T, vs. observed 3-9 keV flux for the same reflection models as 
on the left. 


discussions of this work can be found in Wilms et al. (1998) and Nowak et al. 
(1998b). 


2. X-ray Observations 

We chose to fit the ASCA data with a phenomenological model consisting of a 
multicolor blackbody spectrum plus a broken power law, considered with and 
without a narrow Gaussian line at « 6.4 keV. In all cases, the inclusion of an Fe 
line significantly improved the fit (Wilms et al. 1998). The equivalent width of 
the lines were ~ 40 eV, even for the faintest observation which had 3-9 keV flux 
of ~ 10 -10 ergs s -1 cm -2 . The detectior of an Fe line at such a low flux (along 
with a soft X-ray excess well-modeled by a power law plus multicolor blackbody 
with kT as 150 eV; cf. Wilms et al. 1998 argues against models with extremely 
large coronal radii of 0(1O 4 GM/c 2 ) (e.g., Esin, McClintock, & Narayan 1997). 

We fit ionized ‘reflection models’ after Magdziarz & Zdziarski (1995) to our 
RXTE data. It has been claimed that the best-fit reflection fraction is correlated 
with the photon index of the incident power law, with softer spectra implying 
greater reflection (Ueda, Ebisawa, & Dene 1994; Zdziarski 1998, this volume). 
As shown in Fig. 1, we do not see evident* of such a correlation. A linear fit to the 
correlation is not a statistically significant improvement (Ax 2 < 0.8) over fitting 
with the mean reflection fraction for on;/ of our models. Note that in order to 
account for systematic differences between the PC A (3-30 keV) and the HEXTE 
(17-110 keV) instruments on board RXTE, we fit models to the ins truments 
separately. In joint fits of reflection models, we allowed the normalization and 
slope of the incident power law to be different between the PCA and the HEXTE 
data (cf. Wilms et al. 1998). [The average difference between the PCA and the 
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HEXTE data fits was Ar = 0.12, with the PCA best fit Crab photon index 
being 3 a softer and the HEXTE best fit Crab photon index being 1.4a harder 
than the T = 2.10±0.03 photon index found by Toor & Seward (1974); Wilms et 
al. 1998.] All of our model fits did show, however, that the fainter observations 
were spectrally harder than the brighter observations (cf. Fig. 1.) 


3. Simultaneous Radio Observations 

The first three of our RXTE observations were simultaneous with 843 MHz ob- 
servations taken with the Molongolo Observatory Synthesis Telescope (MOST; 
cf. Hannikainnen et al. 1998), and with 8.3-9. 1 GHz observations taken at the 
Australian Telescope Compact Array (ATCA; cf. Corbel et al. 1998). GX 339—4 
showed a nearly flat (oc j/ 01 - 0 - 2 ) radio spectrum with, typical flux densities of 
7 mJy at 843 MHz. 

An estimate of the minimum size of the radio emitting region can be 
obtained by noting that observationally the brightness temperatures of radio 
sources usually are not larger than 10 12 K, else the electrons will suffer catas- 
trophic inverse Compton losses. The brightness temperature of a uniformly 
bright spherical source is given by (cD / dv) 2 S v /2nk, where d is the diameter of 
the source, D is its distance, v is the observed radio frequency, S„ is the ob- 
served radio flux density, c is the speed of light, and k is the Boltzmann constant. 
GX339-4 is unresolved in the radio, hence a maximum brightness temperature 
places a lower limit on the size of the emitting region. Taking the 7 mJy observed 
at 843 Mz by MOST, we derive 

JS4xl0 12 cm( J ^)«5xl0«je G ( i ^)(^-) . (1) 

This is far larger than the radius of the corona usually assumed in Comptoniza- 
tion models. The derived size for this system is comparable to that inferred for 
the radio emission region of other accreting black hole and neutron star systems 
(cf. Fender & Hendry 1998). 


4. Power Spectral Densities 

All of the timing analysis was performed in an analogous manner to the timing 
analysis that we performed for our RXTE observations of Cygnus X— 1 (Nowak 
et al. 1998a). In Fig. 2 we present the power spectral densities (PSD) in the 
0-21.9 keV band for our faintest and second brightest observation. All but the 
faintest observation showed evidence of an « 0.3 Hz QPO. The characteristic 
‘break frequencies’ of the PSD for the faintest observation were approximately 
a factor of three less than those for the PSD of the brightest observation. The 
faintest observation also showed greater variability amplitude. We found that 
all the PSD could be reasonably well-fit by the sum of a power law plus a 
low-frequency and a high-frequecy broad Lorentzian (cf. Nowak et al. 1998b). 
Below, we discuss the implication of such a fit for phenomenological models of 
the coherence function. 
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Figure 2. Power spectral densities for our faintest (left) and second 
brightest (right) RXTE observation of GX 339-4. Solid line is the best 
fit power law plus two broad Lorentzians. (Fit range is 0.03-10 Hz for 
the left, and 0.1-30 Hz for the right.) Dashed lines show the individual 
components of the fits. 


5. Coherence Function and Time Lags 

We used cross-correlation techniques to compare the variability in the soft and 
hard X-ray bands. We considered lightcurves in four energy bands: 0-3.9 keV, 
3.9-7.5 keV, 7.5-10.8 keV, and 10.8-21.9 keV. The techniques we used are exten- 
sively discussed by Vaughan & Nowak (1997) and Nowak et al. (1998a), and 
references therein. Coherence functions and time lags for the second brightest 
observation are shown in Fig. 3. Note the ‘notch’ in the coherence function at 
~ 2 Hz. This notch may indicate that the decomposition of the PSD into a power 
law plus two Lorentzians may have a physical basis (cf. Nowak et al. 1998b). If 
each of these individual components represents a process that is coherent with 
itself (with constant phase lag between hard and soft variability), but that is 
incoherent with the other two process, then this ‘notch’, as well as other features 
of the coherence can be reproduced, as in Fig. 4 (Nowak et al. 1998b). 

The one clear trend that we have found for the time lags is that the time- 
lag between hard and soft X-ray variabil: ty grows shorter with decreasing X-ray 
flux. We show this trend in Fig. 5. This trend apears to be counter to models 
where the ‘coronal radius’ increases with decreasing flux (e.g., Esin, McClintock, 
& Narayan 1997). Furthermore, considering that the characteristic timescales 
of the PSD are increased (implying lowtr PSD frequencies) for lower flux, this 
trend is counter to energy dependent ‘shot noise models’ (e.g., Poutanen & 
Fabian 1998, this volume). The longer shot timescales required to produce 
lower frequency PSDs should also lead to longer time delays between the hard 
and soft variability. This trend can be countered to some extent, however, by 
also varying the characteristic shot spec ,ra as a function of flux, although it is 
not clear whether the faint and bright sp jctra are sufficiently different to permit 
the required change to the individual shot spectra. 
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Figure 3. Coherence function (left) and time lags (right) of various 
energy channels [squares: 3.9-7.5keV, triangles: 7.5-10.8 keV (coher- 
ence offset by 0.5), diamonds: 10.8-21.9 keV (coherence offset by 1.0)] 
as compared to 0-3.9 keV. Data shown for our second brightest obser- 
vation. Left: Solid lines correspond to unity coherence. Right: Crosses 
are where the soft variability lags the hard variability. Solid line is the 
Poisson noise level for the 10.8-21.9 keV time lags. 



Figure 4. Left: Phase lags between 0-3.9 keV and 10.8-21.9 keV 
variability as a function of Fourier frequency for the faintest obser- 
vation. (All points are hard variability lagging soft variability.) Right: 
[1 — 7 ^ (/)], where 7 ^(/) is the measured, noise subtracted coherence 
function for the same data as on the left. All data has been logarith- 
mically binned over frequencies / — ► 1.4 /. Solid lines are the best fit 
results obtained by assuming that the power law plus two Lorentzians 
fit to the PSD are individually coherent, but are incoherent with one 
another [cf. Nowak et al. 1998b, equations (3) and (4)]. 
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Figure 5. Time lags at Fourier frequencies 0.1 Hz (diamonds) and 
0.9 Hz (squares) as a function of measured 3-9 keV flux. Lags are for 
10.8-21.9 keV compared to 0-3.9 keV variability. 
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